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DEVELOPMENT OF THREE-DIMENSIONAL PRESSURE -DISTRIBUTION 
FUNCTIONS FOR LIFTING SURFACES WITH TRAILING- EDGE CONTROLS 
BASED ON THE INTEGRAL EQUATION FOR SUBSONIC FLOW 

By Alfredo Navarro Crespo* and Herbert J. Cunningham 
Langley Research Center 

SUMMARY 

The problem of calculating subsonic pressure distributions on thin lifting surfaces 
with slightly deflected or deformed control surfaces is studied by the approach of the 
linear-integral-equation procedure. The results are expected to find application to the 
oscillating case. The singularities introduced in the pressure distribution by the deflec- 
tion of a control surface are studied by an analogy between the fluid flow suddenly 
deflected when passing over the hinge of a control surface and the flow also suddenly 
deflected at the sharp (zero-radius) leading edge of a wing with symmetric thickness, 
flying with zero angle of attack at subsonic speed. 

As a result of this study, a suitable set of loading functions (or modes of pressure) 
has been generated to represent the influence of the control surfaces on the pressure 
distribution on the wing. These loading functions can be applied to any wing planform, 
and once they have been integrated through the mentioned integral equation, they produce 
a set of modes of downwash which has a fairly good capability to represent the first-order 
discontinuity that the actual downwash distribution displays along the contours of the con- 
trol surfaces. These resulting modes of downwash should improve the conditioning of 
the matrix equation to which the associated collocation procedure leads, in comparison 
with previously used pressure modes. This effect is very important in the numerical 
calculations. 


INTRODUCTION 

A great deal of work has already been done on the problem of calculating steady or 
unsteady aerodynamic pressure distributions on lifting surfaces in subsonic flow. (See, 
for example, refs. 1 to 20.) As a result of this work, it is possible at present to predict 
with fairly good accuracy the subsonic aerodynamic loads on a thin wing of any planform, 
if the frequency of oscillation of the wing is not too high and the Mach number is not too 
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close to unity. Different improvements of the already well-established calculation pro- 
cedures have been proposed (refs. 12 to 16) and useful work has also been done on the 
case of nonplanar lifting surfaces (e.g., refs. 19 and 20.) 

But there are problems associated with the lifting -surface theory which still remain 
unsolved or are in an early stage for their practical solution. Such is the case of the 
problem of determining pressure distributions on thin wings of arbitrary planform with 
slightly deflected or oscillating control surfaces. An important problem of the latter 
type is the flutter of thin wings with control surfaces. This problem should be amenable 
to linear -theory analysis because just when this instability appears, it is very possible 
that the flow can still be considered as a potential flow and the perturbations as small. 

For many years, methods of calculation based on the lifting-line theory (or the 
strip theory) have been used to determine the aerodynamic loading that the control sur- 
faces produce (refs. 21 and 22). Those procedures, however, have proven to be inade- 
quate for low-aspect-ratio wings and partial-span flaps in common use. 

The reverse-flow theorem (ref. 7) has provided a useful tool for determining gen- 
eralized aerodynamic forces on wings of any planform with control surfaces. (See sec- 
tion 11 of ref. 4.) But its use does not make possible the calculation of the pressure 
distribution itself. 

Many attempts have been made to solve the wing— control -surf ace problem by 
employing the approach successfully used to solve the lifting- surface problem for the 
wing alone. (For the latter, see refs. 1 to 6, for example.) Briefly, the procedure is as 
follows: In the integral equation that relates the usually unknown aerodynamic lifting 
pressure and the usually known downwash distribution, the pressure distribution is 
represented by a truncated series expansion of loading terms (or pressure modes) with 
the terms multiplied by their respective weighting factors, which become the unknowns 
of the problem. Then, the weighting factors are determined by a collocation or a least - 
squares-error procedure. 

In applying this approach to wings with controls the logarithmic pressure singular- 
ity at the hinge line indicated by two-dimensional theory (e.g., ref. 23) has usually been 
accounted for as in reference 24, for example. For partial-span controls, difficulties 
have been encountered in accounting for the ends of the controls and most results have 
not been satisfactory enough to be published. The reason for much of these poor results, 
as experience from the present investigation has shown, is that multiplication of the two- 
dimensional chordwise pressure distribution by separated-variable functions of chord and 
span variables (notably power series) produces individual downwash modes that in linear 
combination have little capability of representing even the statically deflected control. 
Furthermore, the individual higher-order downwash modes tend to be too similar in shape 
and thus contribute to a poor conditioning of the downwash matrix in the collocation 
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process. A recent report (ref. 25) describes results of an effort of this general type. 
Fairly good success was achieved, but some results displayed a noticeably erratic behav- 
ior as the number of pressure modes was increased. 

The conclusion is reached herein that, in general, the two-dimensional chordwise 
distribution of pressure does not reliably represent the pressure distribution for a 
finite-span control. New pressure -distribution functions, developed from three- 
dimensional theory, are needed. 

A new insight to the problem has been given by Landahl in reference 26. The local 
flow at points located near the control surface hinge, side edges, and corners was studied 
by use of the method of matched asymptotic expansions. That study provides very use- 
ful information which can be applied to originate the first-order mode of pressure asso- 
ciated with the control -surface deflection. 

The control surface treated herein is the trailing-edge, flap type with sealed 
leading-edge hinge and no gap at inboard and outboard edges. The terms "control sur- 
face" and "flap" are used interchangeably, Adaptation could be readily made to the 
geometry of leading-edge controls, but the present report treats only trailing-edge 
controls. 

The aim of the present paper has been to originate a set of pressure modes which 
can be applied to any planform with statically deflected and deformed control surfaces 
and, which when transformed through the integral equation previously mentioned, can 
produce a set of downwash modes able to represent adequately, with only a few terms, 
actual downwash distributions and, at the same time, feature strong linear independence 
of the downwash equations. The results are expected to find application to oscillating 
as well as to steady-state problems. 

The way to originate these modes of pressure has been by the use of an analogy 
between the fluid flow suddenly deflected at a control-surface hinge on the upper or lower 
side of a lifting surface and the flow also suddenly deflected by the sharp (zero-radius) 
leading edge of a wing of symmetric thickness flying at zero angle of attack. The down- 
wash distribution on this wing of symmetric thickness is expanded in a series of orthog- 
onal terms. This originates a resulting pressure distribution also expanded in a set of 
terms which can be used to define appropriate modes of pressure needed in the lifting- 
surface case. 

The fundamental or first-order term obtained in this way displays a logarithmic 
singularity which is identical with that found in reference 26. Appropriate modifications 
have been introduced to account for the wing-edge effects as well as the effect of a swept - 
back control -surface hinge, which appears to be very appreciable. The higher order 
modes of pressure have been generated by following a systematic procedure of modifying 
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the first-order term. A parameter E s has been introduced that, if given a suitable 
value, has the effect of improving the downwash distribution produced by the first mode 
of pressure. The suitability of the pressure -mode functions is assessed by plotting 
their associated downwashes and making a visual comparison of these downwashes with 
desired downwashes and with each other, 

Certain auxiliary functions are given in appendixes A and B, and appendix C 
describes the numerical integrations employed. 

SYMBOLS 

A bar over a coordinate variable indicates that it is dimensional, 
weighting factors for wing modes of lifting pressure 

column matrix of subcolumns a nm and b nm 

wing modes (distributions) of lifting pressure (see eqs. (3) and (5)) 
local wing semichord 

weighting factors for control-surface, or flap, modes of lifting pressure 

flap modes (distributions) of lifting pressure (see eq. (3)) 

preliminary exploratory choice for B nm (see eq. (6)) based on two- 
dimensional flow 

second exploratory choice for Bqq (6,v) (eq. (21b)) without the empiri- 
cal modifications that lead to equation (24) 

designation for Bqq(6,tj) of equation (24) when applied to a swept flap, but 
without the sweep effect included in equations (28) 

Bi 0 (M), b’ 2 q(Q,v) incomplete exploratory flap modes of pressure (see eq. (36)) 

c c control-surface chord length 

C nr (M) spanwise symmetric lift distribution functions contained in B nm (0,7?) 

(see eqs. (30) to (35)) 


l nm 


b(7?) 

bnm 
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c H nr (x,y,z), C Tnr (x,y,i) 


dnr 

E S 

F nr( ? ?; x >y) 

F n r (??;x,y) 

Enr(®> 7 ?) 

k 

K 

K 

l 

L A n ( 0 )> l Bj 
m c 

M 

NnrC^j^c? 7 ?) 

n o(e,e c ,v) 

Ap(£,7j) 


perturbation-pressure distribution functions that contain 
and account for leading-edge and trailing-edge singulari- 
ties, respectively, of the thick- wing analogy (see eq. (14) 
and appendix B) 

constant coefficients (see eq. (11)) 

span expansion-contraction factor (see eq. (23d)) 

chordwise integral (eq. (C2)) 


chordwise integral (see eq. (C8)) 
defined with equation (CIO) 

bgw 

reduced frequency, referred to bg, k = — 

kernel function of the integral equation (1), K = K[M,x-£,s(y-??)J 

dimensionless part of K = ■ — — ■ — — as applied to equation (Cl), 

bg s 2 (y -’Vr 

K s K[M,x-|(0),s(y-r7)] 
wing semispan (see sketch 2) 

^(8) chordwise distributions of lifting pressure for wing modes A nm (0,rj) 
(see eq. (5)) and flap modes B nm (9,rf) (see eq. (7)), respectively 

slope of parallel leading and trailing edges of swept control surface whether 
part of a wing or not (see sketch 6 and eqs. (25) and (29)) 

Mach number of unperturbed air stream 

distribution functions in the lift (see eqs. (C5) to (C7)) 

distribution function in the lift (see eqs. (C4)) 

lifting pressure distribution (eqs. (1) and (3)), positive up, 

Ap(£,??) “ P iU,V) ~ P U U,V) 
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p P u (|,??) pressure perturbation on lower and upper sides, respectively, 

of wing surface 

p g (x,y,z) pressure perturbation upon and originated by wing of z-symmetrical 

thickness 

Pjj (x,y,z) local distribution of p g (x,y,z) in vicinity of wing leading edge 
(see eq. (20)) 

P n [ J, P r (^ J Legendre polynomial of first kind of degree n and r, respectively, 
and argument [ ] (see eq. (11)) 

Q r _i( ) polynomial of degree r - 1 and argument ( ) introduced in equation (17) 
and listed in appendix A 

s = l/ b 0 

s e effective value of s (see eq. (23d)), s e = sE s 

S w planform area of full-span wing 

U speed of unperturbed air stream 

w(x,y) downwash distribution at wing surface, positive up (see eq. (1)) 
w 0 (x), w 1 (x)^J 

W 2 (x), W 3 (x),\ illustrative downwash distributions (see sketches 3 and 4) 

w 23 ( x ) J 

w s (x,y,0) downwash distribution on upper side of wing of z-symmetrical thickness 

{wj matrix column of downwash ratios w/U at control points (eqs. (8)) 

W Anm (x,y), w B nm ( x J) distribution of downwash ratio w/U at wing surface 

originated by wing modes and flap modes of pressure, 
respectively 

Wnm( x >y) includes both W Anm (x,y) and W Bnm (x,y) 
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[w] matrix of values of W^ nm (x,y) and W Bnm (x,y) at control points 

(see eqs. (8)) 

W B,lo (X ’ y) ’ W Boo (X > y)> W B00 <X ’ y> 

W B 10 (x ’ y) > W B 20 (X ’ y) 


downwash distributions originated by flap modes of 
1 pressure B B B 0()(M’ 
B' 1O (0,??), and B2 q( 0,??), respectively 


x,y,z dimensional rectangular coordinates (see sketch 1) 

x,y,z dimensionless coordinates, x = x/bq, y = y/l, z = z/bq 

x c (y), x c (y) chordwise coordinate of control -surface hinge line 
x m (y), x ra (y) chordwise coordinate of wing midchord line 


X 


Xn 


distribution function introduced into flap lifting-pressure modes (see 
eqs. (23) and (29)), X = X(6,e c ,i]) 

see equation (23c), X-j = X(tt ,0 c ,v) 


y el , y c j spanwise coordinate of inboard end of partial-span flap (see sketches 2 
and 6) and of left end of full-span flap 


y c 2> y c 2 spanwise coordinate of outboard end of partial-span flap (see sketches 2 
and 6) and of right end of full- span flap 


'cO 


-i 


spanwise midpoint coordinate of flap, y c Q 


y cl + y c2 


1 - M* 


r r (x-l,y,z) distribution function in perturbation pressure (see eqs. (15) and (16)) 


®c = - 


y c2 “ y cl 
5c = y c2' y cl = M 


infinitesimal computational quantity (see eq. (C16)) 
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A sweep angle of wing quarter-chord line, positive for sweepback 

X taper ratio, fy/bp 

X(6 c> r]) see equations (C12) and (€13) 

ju = x T - x c 

M = c c /b Q = p/l b Q 

distribution function (see eq. (17) and appendix A) 

I1q(?]), n(0,6c>n) distribution functions (see eqs. (28)) 
p mass density of unperturbed air stream 

6 angular chordwise coordinate (see eq. (4)) 

6 c (ri) value of 9 at flap hinge line 

6 c q{t]) value of 9 at flap midchord point ^i.e. , at 

6 x {ri) as in reference 1, the value of 9 on span station r] at its intersection 

with the line that passes through the control point (x,y) and is parallel 
to the y-axis 

dimensional and dimensionless {t, = dummy variables for x and x, 

respectively, and with same system of subscripts (see eq. (4)) 

rj,rj dimensional and dimensionless (rj = rj/l) dummy variables for y and y, 

respectively, and with same system of subscripts 

<fi{9,n) see equation (Oil) 

co circular frequency of oscillation 

Subscripts: 

c flap leading edge (hinge line) 
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inboard end of partial-span flap and left end of full-span flap 


outboard or right end of flap 
midpoint or midstation on flap 
wing midchord line 
wing tip 
trailing edge 

wing root or plane of symmetry 

ANALYTICAL DEVELOPMENT 
Formulation of the Problem 

The purpose of this investigation has been to arrive at some practical method of 
calculating the steady pressure distribution on elastic wings with slightly deflected or 
deformed control surfaces in subsonic flow. For this research, the following linearizing 
assumptions were made: a potential flow, an ideal fluid, and small perturbations from 
the main stream. 

These assumptions could seem unrealistic under actual conditions where deflected 
or deformed control surfaces are involved. However, accurate solutions of the three- 
dimensional flow equations, under the mentioned idealized conditions, should be useful 
for the theoretical study and calculation of many aerodynamic and aeroelastic phenomena, 
Such as rolling effectiveness, loss of control effectiveness, hinge moments, and so forth, 
where small deflections of the control surfaces can play an important role. 

In the present approach to the problem, a generalization has been made from the 
common procedure of calculating pressure distributions on lifting surfaces without con- 
trols to the case of control surfaces of immediate interest. That means that deforma- 
tions of the idealized wing surface have been considered, and the development is based on 
the well-known integral equation: 

w fo?y) = M ., H Ap(4,rj) K(M,x-?,s(y-Tjj]d? d ^7 (1) 

u 4?rpU 2JJ S w L J 

where S w is the wing surface. Equation (1) is the steady-state equivalent of equa- 
tion (1) of reference 1. A collocation procedure is used in order to get its numerical 
inversion. 
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The surface integrations have been performed following the technique exposed in 
reference 1 , with proper modifications which take care of the singularities that the con- 
i' trol surfaces introduce in the function Ap (£, 77 ). 

The notations are also analogous to those used 
in reference 1. 

The adopted coordinate system is repre- 
sented in sketch 1 , where S w represents the 
idealized wing surface, which lies very close to 
the xy -plane. The function Ap(|,? 7 ) is defined 
as Ap(£, 7 j) = PjU.Tj) - p u (£,rj), that is, the local 
pressure difference between the lower and upper 
surfaces of the wing. By definition, the function 
is positive when it produces a local lifting force 
directed towards the positive sense of the z-axis. 
Consistent with the adopted coordinate system, w(x,y) is also positive upwards. The 
kernel function K is defined as in reference 1. 



The planform of the wing and control surface considered is represented in sketch 2 
its geometry being fixed by A and the parameters 



x 

Sketch 2 





( 2 ) 


The formulation and the computing program could 
be adapted to a more general planform, in partic- 
ular to tapered controls, but the one programed is 
general enough for the present purpose. 

The main feature that a control surface intro 
duces in the problem of calculating pressure dis- 
tributions on a wing, under the hypothesis of a lin- 
earized flow (M < 1) and idealized conditions of a 
sealed hinge and side edges, is a distribution of 
logarithmic singularity on such a hinge. If this 
singularity is approximated by the usual pressure 
mode series associated with a wing without con- 
trols, slow convergence may result as the number 
of modes is increased, or the series may be found 
to be inapplicable. This circumstance, together 
with other possible problems associated with gaps 
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and possible flow separations, makes it necessary to treat the ease of control surfaces 
with special care. 

Following formulation of the idealized problem through integral equation (1), atten- 
tion focuses on a search for a suitable analytical representation of the function Ap (!,??) 
and on accurate numerical integration over the wing surface S w . For a given downwash 
distribution w(x,y) over the wing, the function Ap(|,rj) should be represented by a 
suitable series expansion which converges as quickly as possible toward the exact solu- 
tion of the integral equation and produces at the same time a good conditioning of the 
matrix of the set of linear algebraic equations to which the collocation procedure leads. 
Achievement of these two properties would mean that the kernel transforms of the pres- 
sure modes would yield a series of downwash modes which should combine to give rapid 
convergence toward a prescribed downwash w(x,y). In order to determine how accurate 
this approximate solution is, the resulting downwash distribution is calculated and visu- 
ally compared with the desired downwash w(x,y). 

In order to represent conveniently the loading due to both wing modes and flap modes 
of pressure, a pressure expansion is chosen of the form 



^ a nm^nm(M) + ^ 


frnmBnm^) 7 ?) 


where the transformation 

h = &&) = (V) - cos 6 


(3) 


(4) 


has been introduced as in reference 1. Although the pressure modes are the products of 

u ■ 

— \Jl - 77^ times A nm (0,r/) and times B nm (0,r/), for convenience in the rest of the 

report they are referred to simply as A nm (0,rf) for wing modes and B nm (0,r/) for flap 
modes and they can be combined additively because the governing equation has been lin- 
earized. The second subscript m is even for spanwise symmetry and odd for spanwise 
antisymmetry of the functions. 


The wing modes can be defined in any conventional form. (See refs. 1 to 5.) As in 
reference 1, the following form was chosen: 

-A-nm (&>V) = 7 ? m I J A n (^ 


where 


La o ( 0) = cot - 
L A n ( 0 ) = 1 


>2n-l 


sin nO 




> (5) 


(n = 1, 2, 


•) 

J 


and even and odd m apply to spanwise symmetry and antisymmetry, respectively. 
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Equations (5) represent the pressure distributions over the complete lifting surface 
without the influence of flap deflection. Since the highest order singularity in the kernel 
function is the same for both steady and unsteady flow, the singularities in the pressure 
distributions can be expected to be the same for both steady and unsteady flow. Hence, 
the pressure modes appropriate to steady flow are expected to be applicable to unsteady 
flow as well. 

The flap modes are developed in the rest of the present report. The coefficients 
a nm an( t b nm are the initially unknown weighting factors in the pressure distribution. 

Preliminary Investigation Based on Separation of Variables 

In preliminary calculations the special case was studied of M = 0 and a flat wing 
of rectangular planform with an aspect ratio of 4 (s - 4, X = 1, and A = 0) and a deflected 
full-span flap hinged at 75 percent of the wing chord. Subsequent discussion in this sec- 
tion related to the inadequacy of the two-dimensional type of loading function is based on 
these calculations. 

In order to gain insight into the problem, begin by assuming the following form for 
the flap loading functions: 

B^(M = r? m L Bn (e) (6) 


(' 


the notation 


*(D, 


B v ‘ l/ ( 0 , 77 ) 
nnr ’ u 


is used to distinguish these as exploratory modes) 


Ij B n (e) = — ^ -(eos 0 - cos 0 c ) n ln 


1 - cos(0 - 0 C ) 
cos 9 - cos 0 C 


(7) 


where 6 C is the coordinate of the flap hinge. For n = 0, 1, and 2, these modes con- 
form to the velocity potentials in reference 23. 

The downwash input w(x,y) considered in the calculations is 0 over the wing 
and 1 over the flap. Eighteen modes of pressure were used to represent the pressure 
distribution: All combinations of n = 0, 1, 2 and m = 0, 2, 4 for both kinds of pres- 
sure modes A nm (0,r? ) and B nm (0,rj). 

The numerical inversion of the integral equation (eq. (1)) relating w(x,y) with 
ApU,? 7 ) was programed on a computer, and the resulting pressure distribution with the 
weighting factors now known was the new input required to recalculate, through equa- 
tion (1), the corresponding downwash distribution. Although this recalculated downwash 
and the initial input w(x,y) coincided at the collocation points as required, the agree- 
ment elsewhere was poor and remained poor through various modifications of the numeri- 
cal chordwise integration procedure used for the modes B ^ (6,r?). 
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This result is mentioned because this kind of representation, inspired like 
B nm (M) ^y the two-dimensional theory, has been for years a common approach that 
has been tried, often without satisfactory results, to analyze the aerodynamic forces on 
three-dimensional (finite-span) control surfaces. 

Neither the use of new flap modes, made with combinations of modes B ^ (0 ,t?), 
nor the use of a least-squares procedure (with 32 control points) for calculating the coef- 
ficients a nm and b nm had an appreciable effect of improving the solution. With the 
least-squares procedure, better shapes of the recalculated downwash distribution were 
obtained in the regions where some number of control points were concentrated, but the 
results remained poor in other regions of the wing. With a limited number of pressure 
modes, a more dense distribution of control points can average on the wing surface the 
good or poor quality of the results but not make them substantially better. Alternatively, 
increasing the number of modes of pressure does not always improve these results but 
can even make them worse if the solution does not converge. 

In the calculations for which the same number of control points and unknown coef- 
ficients a nm and b nm were used, the location of the control points over the wing also 
showed a great influence on the results. An apparently insignificant change in the control 
point locations can produce a drastic change in the shape of the recalculated downwash 
but, in general, without introducing any significant improvement in the total result. 

From these results, equation (7) is concluded to be an unsatisfactory representa- 
tion of the pressure distribution. 

Another possible source of trouble which can considerably influence the results 
is inaccuracy in the numerical integrations. Accurate techniques can be used to perform 
the numerical integration of functions which contain a logarithmic singularity, such as 
the Berthod-Zaborowski quadrature (a weighted Gaussian type of quadrature, see, e.g., 
table 9 in ref. 27). But in the present case, the kernel function introduces an additional 
problem (it approaches a chordwise step function in the vicinity of the collocation point; 
i.e., as y - 77 — 0 ) which is difficult to treat and can be a source of important errors, 
mainly when the control point is located very near the control-surface hinge. The numer 
ical surface integration of a logarithmic singularity which terminates within the region 
of integration (such as that at the hinge of the usual control surface with a span smaller 
than the wing span) also presents some special difficulties. This is the kind of logarith- 
mic singularity that the functions B nm (0,r?) possess, as will be shown in a subsequent 
section herein. Nevertheless, experience gained during the preliminary calculations 
has given confidence in the accuracy of the integration procedure adopted (also described 
below) even though the quality of the pressure representation was not satisfactory. 

With the purpose of gaining better insight into the reasons why this representation 
of the pressure distribution was so poor, the downwash distribution which corresponds to 
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each mode of pressure used in equation (6) was studied. These distributions are called 
wing modes of downwash WA nm (x,y) and flap modes of downwash w B nm ( x >y)- 

Each downwash mode WA nm (x,y) or Wg nrT1 (x,y) comes from its related pres- 
sure mode A nm (d,r]) or B nm (0,?7), respectively, through the basic integral equation. 

The function becomes represented by the expansion 

j a nm w An m ( x >y) + ^ b nmW Bnm ( x ,y) (8a) 

n,m n,m 

where n = 0, 1, . . . n max and m = 0, 1, . ... . m max . (The maximum values n max 
and m max can be different in each summation.) 

The evaluation of this truncated expansion at some particular points (control points) 
on the wing surface gives the matrix expression 



from which it is possible to calculate the values of the weighting factors a nm and b nm . 
Good convergence of this downwash expansion and good conditioning of the matrix [w] 
are qualities which are sought when a collocation procedure is adopted. 

A convenient and possibly ideal choice of pressure modes A nm (6,rj) and Bnm(S,V) 
to be used in the pressure expansion would be that which originates a set of downwash 
modes W nm (x,y) that not only is orthogonal but also converges rapidly (with a few of 
the lowest order modes) to a good approximation of the downwash for actual deflecting 
and deforming controls. Furthermore, the higher order terms should produce small and 
decreasing corrections to adjust closer to the desired downwash the approximation that 
results from the truncated expansion. 

These qualities are particularly important in cases where control surfaces are 
involved. The deflection of a control surface introduces a first-order discontinuity in 
the function w(x,y); that is, the downwash is a nonzero constant on the control surface 
and zero on the rest of the wing. This elementary mode of downwash accounts for the 
deflection of a control surface in the function w(x,y). Furthermore, for the case in 
which the control surface behaves as a rigid structure , this downwash mode is the only 
one which introduces the logarithmic singularities in the pressure distribution. Sketch 3 
illustrates this first flap mode of downwash. 
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w(x) = w Q (x) + w x (x) - w 23 (x) 

Sketch 3 

The sketch represents a cambered wing section, taken in a region where a control 
surface has been deflected, and the chordwise downwash distribution which corresponds 
to this section. Also shown is the separation of the actual downwash w(x) into three 
modes of downwash: Wq(x), w^(x), and w 23 (x). 

The mode w^(x) is the first flap mode of downwash. The mode w 23 (x) can be 
further separated into two more elementary modes, as is illustrated in sketch 4. 



Sketch 4 
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( 9 ) 


Thus, the downwash w(x) can be separated into the four elementary modes 

w(x) = Wq(x) + Wj(x) - Wg(x) - Wg(x) 

The mode Wg(x) is one of the kind called wing modes of downwash WA nm (x,y). 

(it could be any chordwise third-order mode, represented by the general expression 
WA 3m ( x ,y)-) The modes Wj(x), w«(x), and Wg(x) are flap modes of downwash, which 
are defined by W Bnm (x,y). (They could be any one of WB 0m (x,y), W Blm (x,y), or 
W B2m( x >y)> respectively.) The mode w^(x) represents the first mode W B Qg(x,y), 
when the downwash distribution on the control surface is a constant in the spanwise sense. 
This mode, together with a smaller contribution of the mode WBg^(x,y), which becomes 
important in cases where the control surface can experience an appreciable torsion, may 
be the only component of the downwash distribution w(x,y) responsible for the logarith- 
mic singularity of the pressure function in many practical calculations. 

Clearly, it would be useful to develop in the pressure expansion a term Bqq(6,t]) 
able to originate a downwash distribution WgQg(x.,y) as close as possible to this ideal 
mode w^(x). If none of the modes W Bnm (x,y) appear to be a good approximation to this 
ideal shape, a high number of modes may be required to represent adequately the actual 
downwash and, consequently, the corresponding pressure distribution. This characteris- 
tic will also aggravate the additional problems of convergence and good conditioning of the 
martix [wj . 

The "optimal location" of control points actually in use (see refs. 2 to 6) is based 
on a separation of variables which does not always yield significant improvements in 
practice in comparison with arbitrary location, at least for low frequencies and simple 
downwash modes (see refs. 17 and 18, for example), and which does not apply to the case 
of control surfaces. Anyway, the results of lifting-surface calculations depend more 
directly on the capability of the downwash truncated expansion to approximate the func- 
tion w(x,y) than on the optimizing of the control-point locations. That does not neces- 
sarily mean that the search for a realistic optimal conditioning of the matrix [wj is 
unimportant for improving the results in calculations of wings with or without control 
surfaces. But it must be kept in mind that this optimal location depends essentially on 
the actual shapes of the modes W nm (x,y), its validity being dependent on the capability 
of these modes to represent the actual downwash distribution on the wing. Furthermore, 
for the analytic expressions and integration procedures used herewith, control-point 
locations must be selected that are not too close to regions where the modes W nm (x,y) 
can suffer appreciable distortions due to numerical errors; such regions are wing leading 
and trailing edges and control -surface boundaries. It is also recommended to avoid, as 
much as possible, the nodal lines of the modes W nm (x,y) and to spread the control points 
on the regions where these downwash modes are widely different among themselves. 
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In figure 1(a) is plotted the first downwash mode w|j^(x,y) obtained by the early 

calculations mentioned for a full-span flap on the wing with the parameters listed. (The 
superscript (1) is used to distinguish these modes as exploratory.) Figures 1(b) and 
1(c) show two higher order modes W^^(x,y) and W^^(x,y). The mode Wg^(x,y) 

appears to be an unsatisfactory approximation to the ideal shape w^(x) in sketch 3, 
and the higher modes can do very little to improve this approximation. The modes 
W B2o^ X ’ y ) anc * Wg^(x,y) are too much alike for numerical purposes, and this is a 

circumstance which affects adversely the conditioning of the matrix [w] . 

Figures 1(a), 1(b), and 1(c) also show for comparison the wing downwash modes 
WAoo^y)’ w A 2 o( x >y)’ anc * W A30 ( x ,y). The contrast between the good conditioning of 
these wing modes and the poor conditioning of the flap modes is evident. This is the rea- 
son why the results of calculations of pressure distributions on wings without control sur- 
faces are much better than those on wings with control surface. 

These particular results show that any combination of wing and flap modes shown 
would have a poor capability for representing the desired downwash w(x,y) for a 
deflected flap and would cause a worsening condition of the matrix [w] if an attempt 
were made to improve the results with additional higher modes from equation (7). 

In contrast, the two-dimensional wing theory from which the pressure modes 
A nm (6,ii) and B^( 0 ,i 7 ) are selected provides a good basis for comparison of these 
results. It is known from the existing literature (see, e.g., ref. 23) that for the analo- 
gous case of a thin planar airfoil with a deflected flap, a combination of only the two first 
modes Laq(^) and LgQ(0) is enough to represent the exact solution of the lifting pres- 
sure in steady flow. The term L^q( 0) originates a downwash mode which is a constant 
on the total airfoil, and the term Lbq( 0) originates a downwash which takes a constant 
value on the flap region and another substantially different constant value on the rest of 
the airfoil. That should be an ideal goal for the three-dimensional case too. But figure 1 
shows that this goal has not been achieved in the early calculations based on equation (7). 

These results make evident the need for new pressure modes Bnm(0,7) able to 
originate a set of better conditioned downwash modes WB nm (x,y). 

Search for Surface Loading Functions 

Because of the desired nature of the downwash mode WgQQ(x,y) (as close as pos- 
sible to a constant value on the control surface and a different constant value on the rest 
of the wing) , one cannot expect to get a general orthogonal set of modes W nm (x,y) that 
includes both wing and flap downwash modes. But it appears to be possible to get wing 
modes which are already fairly close to a set of orthogonal functions (by combining the 


17 



chordwise pressure distributions L^ n (6) with proper spanwise distributions), then to 
find flap modes WB nm (x,y) which also are nearly an orthogonal set. The problem is 
how to find a set of pressure modes able to originate these downwash modes. 

The following section presents the development of pressure modes based on an 
analogy between (1) the flow perturbations over a nonlifting wedge-shaped wing with sym- 
metric thickness distribution and (2) the perturbations over a lifting deflected flap. 


U 


*cl 


y.c2 




Sw 


Wing section 


Analogy of nonlifting wing with thickness .- It is hypothesized here that for a non- 
lifting finite -planf or m wing with symmetric (with respect to the xy-plane) nonzero thick- 
ness distribution and a sharp leading edge, those fluid particles traveling infinitesimally 

close to the xy-plane and suddenly 
deflected in their trajectories when 
passing by the leading edge are subjected 
to physical conditions very similar to 
those particles traveling close to the 
upper side of a lifting surface and being 
suddenly deflected in passing the con- 
cave corner at the hinge of an upward 
deflected control surface. From linear 
theory, the effect of a convex corner 
would be equal and opposite in sign. 

Such a nonlifting wing with symmetric 
thickness and sharp leading edge is 
illustrated in sketch 5. On the basis of 
this hypothesis, it is expected that the 
solution for the flow past the leading 
edge, including the corners, of the wing 

of sketch 5 can provide useful information on the analytical form of the local pressure 
distribution in the neighborhood of a control-surface hinge and edges. 

If w s (x,y,0) is the distribution of simple sources that matches the downwash dis- 
tribution on the upper surface of this wing (sketch 5) and p g (x,y,z) is the associated 
pressure perturbation, then for the steady wing they are related, according to linear 
theory, by the expression 




Wing planf orm 


Sketch 5 


p s (x,y,z) * - jl rr 

s 2 tt ax 


w s (i,b,o) 


V(x - l) 2 + 0 2 [(y - ri) 2 + z 2 ] 


d| drj 


( 10 ) 


The advantage of this approach becomes apparent. Only simple sources need be con- 
sidered, and the desired pressure is directly obtainable rather than being contained within 
an integral of an integral equation. 
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Consider a function w s (i,?j,0) defined on the wing surface by the series expansion 
of orthogonal terms 


A 0) = u £ d nr p n r|(5 c0 -{)lp r [6 c (y c0 ^) 

n,r L J 


where the quantities d nr are constant factors whose numerical determination is of no 
concern here, and 


(*02 > *cl) 


p = x T - x c 


y c2 “ y cl 


X T + x c 
x c0 = 2 


- y C 2 + y d 

y c0 2 


The functions P n and P r are Legendre polynomials of arguments ^ x c q~IJ anc * 


^ c ( y c0"^)’ res P ectivel y 5 an d can be defined by 




r\/2r-2j' 


3 A *• 


(12a) 


(where [r/2] indicates the integer part of r/2 and f is a dummy argument), or 
alternatively by Rodrigues’ form 


^rW = ^0 2 - l)r ] 

This downwash (eq. (11)) originates a pressure 
P s g,y,i) - - ^ / dnr £ T P„ft(i c0 -4)]di ^ t/ 2 


■[ § c(y C Q-^] 


y cl \(x - i) 2 + /3 2 [(y - V) 2 + i 2 ] 


The integration with respect to rj can be carried out in closed form so that the pressure 

is then represented by the expansion 

m 

p s (5,y,z) = - ^ £ d nr rC Tnr (x,y,i) - C Hnr (x,y,z)j (14) 

n m V* • 
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where Ct and Cu are the values that the indefinite integral 
A nr n nr ■ 



9 /- 

J p " 

p( X cO^)_ 


£ r r (x-I,y,z)d! 
ax 

takes at the trailing edge | = x^, and leading edge \ = x c , respectively. 
The function T r (x-|,y,z) is defined as 


( 15 ) 


- - -x -i ]l& ■ + p 2 (y - y c i) 2 + 52 - p(f ■ y c i) 

r. x-f v 7 - 1 In J Cl/ J V Cl/ 

(16) 

(T 7 Q 

'<* - 1) 2 + S 2 [(y - y c2 ) 2 + i 2 ] - 0 (y - y c2 ) 

When r s 1, 



r r (x-I,y,z) = \j& - 1) 2 + P 2 

_(y-yci) 2 + i2 ] , 3r-i(y-yc!) 



- J(x - |) 2 + /3 2 (y - y c2 ) 2 + z 2 Q r _ 1 (y - y c2 )+ v r (l)r 0 (x~S,y,z) (17) 


where Q r _j is a polynomial (not Legendre's second kind) of degree r - 1 and argu- 
ment y - y cl or y - y c 2 , and v (|) is a function of | (also of x, y, and z) and 
does not depend on r). Both can be calculated by the method of undetermined coefficients 
through the identity 


[ 5 c(y c o-^)] s - i) 2 + p 2 (y - ^) 2 + z^jQr-i(y-^) + P 2 (y - n)Q r _i(y-t?) + v r (l) 


by matching on left and right sides the coefficients of each power of rj, the coefficients 
themselves being functions of x, y, and z. In equation (18), 

%.^-ri) = Q r _ 1 (y-’7) (19) 

Expressions for Q r _j(y-77) and y (1) for r up to 5 are given in appendix A. 

The local distribution in the neighborhood of the wing leading edge (which simulates 
a control-surface hinge) is given by the expression 


* p Hs (x,y,z) = ^ Yj d nrCHnr^’y^) ( 20 ) 

n,r 

By taking the limit z — 0, the local pressure distribution Pjj g (x,y,0) in the xy-plane is 
obtained. This function also represents the local pressure distribution which appears on 
the upper side of a flat lifting surface in the neighborhood of the hinge of an upward 
deflected control surface, when the downwash is zero on the wing and w s (x,y,0) on the 
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control -surface region. The terms Cjj nr (x,y,0) of the resulting expansion are then 
pressure distributions associated with an orthogonal set of downwash modes. They con- 
tain the appropriate singularities that the hinge of a control surface introduces in the 
pressure distribution, and with the usual corrections to account for the effect of wing tips 
and leading and trailing edges (as in ref. 1), they can serve to originate the flap modes of 
pressure B n m (6,V) that are sought for the lifting-surface case. 

The terms C H()0 (x,y,0), C Hl0 (x,y,0), C H20 (S,y,0), C Ho i(x,y,0), and CH 02 (x,y,0) 
are given in appendix B. Since the main effect of deflection of the control surface is given 
by the mode Bqo (Q,v) and since the higher order modes B nm (0,?7) serve mainly as cor- 
rections of the distortion that the wing modes A nm (0,??) introduce in the downwash dis- 
tribution on the control-surface region, only a few flap modes of pressure B nm (0,?7) are 
expected to be needed in most of the practical calculations. The relationship between the 
indices m and r is given subsequently in the development of higher order pressure 
modes, where it becomes necessary for m > 0. 


The term CjjQQ(x,y,0) has the expression (in dimensional coordinates) 


C Hoo (x,y,0) = I in ) M tiy 

|/(x - x c ) 1 2 + /3 2 (y - y c2 ) 2 - /?(y - y c2 ) 


(21a) 


which except for notation is the same as the local solution found by Landahl (eqs. (33) 
to (35) of ref. 26, with reduced frequency set equal to zero). In reference 26, the pres- 
sure distribution for an oscillating planar rectangular control surface was developed by 
the method of matched asymptotic expansions. 

One way of introducing the correct wing-edge behavior is by multiplying 

CHoo(x,y,°) by the factor 

I (x - x L )(x T - i) 

f(*c ~ %)( x t " x e) 

as in reference 26. The resulting flap pressure mode is denoted for a wing with partial- 
span flaps (as in fig. 6) by 


1 sin D j cos 6 - cos g c) 2 + ^ s2 (l ^ I - ?cl) 2 - Ml ” 1 - ?cl) 
y(cos 6 - cos e c y + ^ 2 s 2 (| r] | - y c2 ) 2 - 0s(j?7 j - y c2 ) 

where the variables x and y have been changed into 6 and 77, respectively, and a 
superscript (2) is used to differentiate this mode from the finally adopted Bqo( #,??). 
For a full-span flap (as in fig. 2), y - = -1, y c2 = 1, and J 77 1 is replaced by 77. 
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In figure 2, is shown the downwash distribution (x,y) that the mode B^(9,rj) 

originates when applied to a square planform wing with a full-span flap rj c ^ ~ -7 7 c y = 1.0. 
A good characteristic of this downwash distribution is that it has at the hinge nearly a 
first-order discontinuity in the chordwise sense, and this discontinuity has a nearly con- 
stant magnitude along the span. This quality is found in all the numerical calculations 
performed with the mode That means a good approximation to the ideal first- 

order discontinuity which is sought in the mode Wg (x.,y). But the mode w bqq( x >Y) 

also presents, for a reasonable amplitude of the discontinuity, the undesirable qualities 
of a high slope chordwise and a too high curvature in the spanwise sense. These qualities 
produce a general shape which is not close to the flat distribution (with a first-order dis- 
continuity on the control-surface contour) desired for the mode Wg^(x,y). 

There is another difficulty involved in the mode Bqq (0,77): If s increases in such 
a way that the wing planform tends to become a two-dimensional wing with a deflected 
flap, the values of the pressure distribution on the wing also increase and tend to infinity 
in the limit. This result does not make physical sense and indicates that the analytical 
expression of the mode must be corrected. 

This difficulty would not be present if, together with each term Cji nr (x,y,0) in the 

expansion of the local pressure distribution p™ (x,y,0), there had been included the effect 

n s 

that the corresponding term CT nr (x,y,0), which also appears in the total solution 
P s (x,y,0), introduces in the limit when the magnitudes y^ and y c 2 increase in the 
mentioned way. 

Empirical modifications of pressure modes .- An attempt was made to improve the 

mode B qq (#,??), and all the following higher order modes, by introducing the following 
empirical modifications to equations (21): 

(1) For use on the lifting case, each term Cn nr (x,y,0) has been "corrected" with 
a function which behaves as C'i' nr (x,y,0) (from eq. (15)) when y c j — -°° and y C 2 — 00 
but which obeys the trailing-edge Kutta condition rather than having the logarithmic 
singularity that each term C'T nr (x,y,0) has on the trailing edge for the nonlifting wing 
of symmetrical thickness. 

(2) The variable x always appears in the combination x - x c in each term 

CH nr ( x >y>°)- 

Introduction of the same transformation as in equation (4) leads to 

£ - Z C (V) = - £—( c °s 9 - cos <9 C ) (22a) 
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where 


Note that in the vicinity 



9 ~ 9 C the quantity 


cos 6 - cos 9c | ~ 6 - 6 C sin 6 C 


(22b) 


In order to cause the functions B nm (0,?]) to approach the two-dimensional pressure 
distribution of equation (7) (for nm = 00) as the wing planform tends toward a two- 
dimensional wing plus flap, wherever | - £ c (?]) appears, it is replaced by 


-x(M c,v) 


b_ 

b 0 


l - cos(e - 9 C ) 
cos 9 - cos 0 C | 


2 sin^g,. 


(23a) 


which is asymptotic to | - £ c (v) as £ — £ c (??); that is, in the vicinity 9 ~ 9 C 

- i*— — — TT 2 sin2 0 c ~ i 0 “ 0 c| sin e c (23b) 

jcos 6 - cos 0 C | c * 

which is the same asymptotic behavior as that expressed by equation (22b). In addition 
the following definition is introduced: 


X T = X(tt,9 c ,v) = - ^ 2 sin 2 0 e (23c) 

(3) The semispan-semichord ratio s is replaced by an effective value 

s e = sE s (23d) 

where E s is a span expansion-contraction factor that can be a constant or a span vari- 
able. The function of E s is to make relatively stronger (Eg < 1.0) or weaker (Es > 1.0) 
the induced effect of the control-surface tips. A value of E s greater than 1.0 has the 
effect of putting the control tips farther away from the rest of the control surface and 
wing. The effect of the control tip is thus attenuated and the pressure distribution in the 
central region of the control surface is closer to the two-dimensional distribution. When 
E s is given an appropriate value, it has a very noticeable effect of improving the down- 
wash distribution Wg^(x,y). This value depends mainly on the wing and control-surface 

configuration, but it is also affected by the Mach number and possibly by the reduced fre- 
quency considered in the analysis of any unsteady case. The influence of this parameter 
and its order of magnitude for several configurations, as well as the influence of the Mach 
number, will be discussed further in the section entitled "Results and Discussion." Since 
no practical relation for predetermining the appropriate value of E s is yet known, some 
numerical exploration in each particular case is recommended. The first mode Bqo (G,rj) 
can be programed with E s as an input parameter, and then several values of E s can 
be tested to pick the one which gives the best downwash mode Wbqq(x,y). 
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With these three empirical modifications to equation (21b), the first flap mode of 
lifting pressure for an unswept flap becomes 


BooCM) 


fe T 2 + jS 2 s e 2 | 

'jc2 - l”l 

f + 0s e (y C 2 - M) 

]/x T 2 + j3 2 S e 2 

(y c i - Kf + fsMjd - hi) 


ijx\e,e c , v ) 

+ £ 2 s e 2 (y c i - M) 2 + £ s e(y c i - M) 

|/x 2 (0,e c ,r?) 

+ /3 2 s e 2 (y c2 ->[ 

f * P s e(y c 2 - M) 


(24) 


where the absolute value | r? | is taken in order to make the expression applicable to both 
right- and left-wing semispans, and where the first factor in the argument of the log func- 
tion (the one containing X-p) is the correction factor resulting from modification (1) 
above. The multiplication of the log function by sin 0/sin 0 C as in equation (21b) is not 
needed because with the given modifications the function Boo of equation (24) has the 
desired zero lift at wing leading and trailing edges. Also note that 

lim Bqo(M) = 4'ln X(e,6 c ,v) 
s e y c i— 00 p 

s ey c 2-°° 

which represents the two-dimensional pressure distribution. (See eqs. (6) and (7).) 

Pressure modes for swept controls .- As will be described in the section "Results 

is produced a downwash mode WgQQ(x,y) 
fairly close to the ideal function that is desired. 
But this pressure mode has been deduced under 
the assumption of a control-surface hinge 
placed perpendicular to the direction of the 
main flow. The numerical calculations have 
shown that when the hinge sweepback angle is 
large, the downwash distribution WBoo( x >y) 
becomes very distorted if the sweep is not 
accounted for, and a more general mode of 
pressure is needed that does account for the 
sweep. 

This more general expression has been 
obtained by returning to the analysis of a thin 
wing of symmetrical thickness and a sweptback 
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Sketch 6 
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leading edge, as represented in sketch 6. The leading edge is represented by the straight 
line 

x c (y) = x c (0) + m c y (25) 

of slope m c . 

A constant downwash distribution has been considered, and the more general func- 
tion CHoQ(x,y,0) takes the form analogous to that for the unswept surface (eq. (21a)): 


c Hoo( 5 >y> 0 > = 


^3 2 + m c 2 


f j2 + m c 2 ij[x - 5 c(y cl ) 

] 2 + P 2 (y - y cl ) 2 - 

(s 2 * »C 2 ) 

(y - y cl ) - m c[x - x c (yjj 

\Ip 2 + m c 2 £x - z c (y c2 ) 

2 - ^ 2 (y - y C 2) 2 - 1 

{? * me 2 ) 

(y - y C 2 ) - m e[x - s c(yj] 


(26) 


From the relations 

x - x c (y c ]) = [x - x c (y)] + |^ c (y) 

x " 5 c (7c 2) = [* ■ ^c(y)] + |~x c (y) 
x c (y) - ^c(y c i) = m c(y - y c i) 
x c (y) - x c (y c2 ) = “ C (y - y c2 ) 


- *cfoi)] 

- M 5 c2)] I 


J 


(27) 


and the use of the same modifications that were made to equation (21a) the definitive func- 
tion Bqo(M) to be used in the definition of the flap modes of lifting pressure is 
obtained: 

B 00 (M) = , 1 : lnfn o (77)n(0,0 c ,?7)l (28a) 

Y/3 2 + m c 2 


where 


||[X T + m c s e(y c 2 ‘ 1 7 1)] 2 + ^ s e 2 (y c 2 "hi) 2 + s e^ 2 + m c 2 (y c2 “ M) 

n 0^) = -- — — 

j|[ x T + m c s e(y c l - h \f + ^ s e 2 (y c l ■ hi) 2 + s e |//3 2 +m c 2 (y cl - | r?|) 


. m c v 

+ An 




•2 * m _2 


m r 


fi 


•Xn 


|2 + m 2 


(28b) 
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(28c) 


Ti(e,e c ,v) = 




|x(Mc>n) + ni c s e (y c 2 - hi)] 2 + ^ 2 s e 2 (y C 2 - H ) + s eP + m c 2 (y c2 - 


7? + 


m c 


+ m c 2 ) 


L(8,0c>ri) 


and the log function approaches zero as 6 — 0 or ir at wing leading and trailing 
edges. As s e becomes indefinitely large, Bqq( 6,77) approaches the two-dimensional 
distribution. 

For the particular kind of planform that is considered, that is, with the flap hinge 
line parallel to the trailing edge, 

X(M c ,v) = [l - h|(l - 2 sin 2 6c] 


6 C = 7r - arc cos 


and 


1 - 


cos 6 - cos 6 C | 
1 




i - hid - x) 




m c = tan A - —(1 - A) 

iiS 


(29) 




When the slope m c = 0, equation (28a) reduces to equation (24) for the unswept wing. 

On the control-surface hinge, the function Boo(8,ri) has a logarithmic singularity 
which disappears beyond the tips of such a hinge. It also has the proper leading- and 
trailing-edge behavior, and is an analytic function on the rest of the wing surface. 

Higher order pressure modes .- In the development of the higher order pressure 
modes, two circumstances must be accommodated. One is that there is an existing usage 
of the second subscript m of B nm (6,r?) as even for spanwise symmetry and odd for 
spanwise antisymmetry. The other is that the subscript r on the spanwise varying 
Legendre polynomial ?r j^c(y c o~^)] introduced into equation (11) can take successive 
values r = 0, 1, 2, . . ., and so forth. These circumstances are accommodated by the 
following definition: 

Bnm(M) = t? C n r(M (r = 0, 1, 2, . . .) (30) 

where p = 0 for spanwise symmetry, p = 1 for spanwise antisymmetry, and 

m = 2r + p (31) 

The higher order functions Cnr(0,?7 ) which appear in the definition of the flap 
modes of pressure B nm (6,r]) can be obtained by introducing the above-mentioned 



modifications in the higher order terms Cji nr (x,y,0) listed in appendix B. Then, since 

Bqo(M) = Coo (M)> 


3io<M) = C 10 (M) = |(y c2 - 1’']) 1 " 


\jx 2 (e,e c ,n) + /3 2 s e 2 (y c2 - |ij|) 2 - x(e,e c ,v) 


fr * S 2 s e 2 (y c2 - 1»|) 2 - x T 


|(yci-i ,) l) ln 


p 2 (e,e c ,T|)+ )3 2 6 e 2 (y el - |i)|f - x(s,8 c ,t,) 


v. COS 0 - COS 0„n 

+ 55 H Coo(M) 


(32a) 


where 


and 


0 cO = arc cos(^ m (,) - £ c0 0?)j 


^c0 (?7) = 


Z T (v) + S C <T7) 


(32b) 


For the swept tapered wing with untapered flap which has been considered, 0 q becomes 


0 cO = ff - arc cos /I - 


M 


Similarly, 


\_ 2[l-|rj|(l-X)] 


(32c) 


S / h \2 1 - cos(e - 0 q) . + (Ts e “ij , - mo - 

. c !0( ,„ . ^ =7^- • I-D- 


.^y % _ |„|) 1 „ l/ x2 ( g ’ g cm) + A e 2 (y c i - hi) 2 - x(6,9 e ,7 ? ) 

^x T 2 + /3 2 s e 2 (y cl - hi) 2 - x T 


rb ' 2 /c Qs e - cos e cQ \ 2 
v b oj \ . p " / ’ 12 


(33) 
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|/x 2 (e,0 c ,7?) + /3 2 s e 2 (y c2 - hi) 2 - |/x T 2 + /3 2 s e 2 (y c2 - | r^) 2 


- + ^ 2 s e 2 (y cl - hO 2 - ifc 2 + ^8^(y c i - M) 2 

+ 166 c (y c0 -|r ? |)c 0 o(M) (34a) 


where 


and 


y c o 


6 e = 


y C 2 + y e i 




y C 2 - y c i 


(34b) 


B o4 (0,r?) = C O2 (0,?7) 


100 


36, 


-(y c o - hi ) 2 


+ 75 


36 r 


4/3 2 s 2 


X 2 (M c ^) 


C 00 (M (35) 


These five modes may be enough to represent the influence of the control surfaces 
on the pressure distribution in many practical calculations. But if still higher order 
modes are needed, it would be possible to generate them by calculating the appropriate 
functions Cfj nr (x,y,0) in the way which has been indicated and applying the mentioned 
modifications to translate these functions to the lifting pressure distributions C nr (0,rj). 

The functions Cjj nr (x,y,0), which come from calculations in which sweepback is 
not considered, are used to generate the general higher order modes. This way is chosen 
in order to simplify somewhat the analytical expressions obtained for the higher order 
functions C n r(0,ij). Anyway, an important influence of the sweepback effect is introduced 
in all these functions through the first-order mode Cqq(6,t}). 


RESULTS AND DISCUSSION 


The results presented in figures 1(a), 1(b), and 1(c) were discussed in the section 
"Preliminary Investigation Based on Separation of Variables." The rest of the results 
are now presented and discussed. Figure 2 applies to a square planform wing (aspect 
ratio of 1.0) with full-span flap hinged at 75 percent of the chord. Chordwise distribu- 
tions of downwash are plotted for the three span stations indicated. At the top is the 
downwash originated by the pressure mode in equation (21b), and at 
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the bottom is from pressure mode Bqo(0,t 7) in equation (24) with E s = 3/2. 

Both distributions display the desired characteristic of nearly a first-order discontinuity 
(finite jump) chordwise at the hinge line and the magnitude of the jump is nearly uniform 
spanwise. Of the two distributions, Wjjqq has the slight advantages of lower slope 
chordwise (away from the hinge) and less variation spanwise. 

Figure 3 applies to the same configuration as figure 2 and the effects are shown of 
decreasing the factor E s (eqs. (23d) and (24)) to three values successively lower than 
that at the bottom of figure 2 . As E s is decreased, at first the chordwise slope of the 
downwash decreases both forward and aft of the hinge, but for the lowest value of E s 
(3/20), some irregularity appears just forward of and aft of the hinge. For this range of 
results, an intermediate value of E s most nearly gives the desired uniform step-function 
discontinuity from wing to flap. 

Figures 4(a) to 4(e) apply to a spanwise symmetric rectangular planform wing (the 
half-span is sketched on the figure) with an aspect ratio of 4.0, partial-span controls from 
0.3 to 0.7 of wing semispan and hinge at wing 3/4 chord. Chordwise distributions of 
downwash are shown at each of three semispan stations. In figures 4(a) to 4(e), E s has 
the values 3/2, 3/4, 1/2, 3/8, and 1/4, respectively. The figures show how effective the 
parameter E s can be; figures 4(b) and 4(c) show that E s in the neighborhood of 3/4 
to 1/2 produces a downwash distribution fairly close to the fundamental (first) flap down- 
wash mode that is desired. Figure 5 shows a thr.ee -dimensional view of the downwash 
distribution for E s = 1/2. 

Figure 6 shows for comparison the downwash distributions W^^(x,y) originated 

by the pressure mode B^(e,rj) of equation (21b) for the same configuration as in fig- 
ure 4. The positive chordwise slope displayed is associated with the multiplying factor 

sin 6 in B^(0,?]). 

Figures 7(a) and 7(b) show for the configuration in figure 6 the effect on W^qq of 
increasing the Mach number from 0 to 0.5 and 0.8, respectively. The value of E s was 
increased to 3/5 for what is judged a better quality of the downwash distribution than that 
in figure 6. For M = 0.8 the distribution is higher at the outboard station y = 0.8 than 
for M = 0.5, and a further variation of E s might have disclosed a better overall choice 
than the value 3/5. 

Figures 8(a), 8(b), and 8(c) apply to an untapered wing with a flap like that of fig- 
ures 4 to 7 except that the wing is sweptback to 45° by shearing (full-span aspect ratio 
remains 4). Downwash distributions are plotted against local chordwise coordinates for 
the same three semispan stations. In figure 8(a), the downwash was obtained by 

applying the pressure distribution Bqo of equation (24) that contains no provision for 
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sweep. The value of E s is 1/2. The need for improvement is apparent. In fig- 
ure 8(b), the results of using Boo °f equation (28), that does account for sweep, are 
shown for E s = 1/2. The effect of increasing E s further to 3/5 is shown in fig- 
ure 8(c), where the downwash is more uniform at y = 0.8 than in figure 8(a). 

Figures 9(a), 9(b), and 9(c) apply to the same parameters as 8(a), 8(b), and 8(c), 
respectively, except that Mach number has been increased from 0 to 0.8, and in fig- 
ure 9(c), E s is 0.556. In comparison with the results for M = 0, the downwash dis- 
tributions for M = 0.8 are noticeably worsened but are still fairly good in figure 9(c). 

Figure 10 shows downwash distributions WbqqCXjY) f° r a rectangular wing with 
an aspect ratio of 2, for M = 0 and E s = 3/5, that can be compared with the distribu- 
tions for the wing with a higher aspect ratio in figures 4(b), 4(c), and 5. A strong simi- 
larity is apparent. 

Figures 11(a) and 11(b) apply to a tapered, swept wing with partial-span constant- 
chord flaps. The taper ratio X is 0.6, the quarterchord sweep angle A is 35°, the 
ratio s is 2, the ratio p of flap chord to root half-chord is 0.4, and the parameter 
E s = 0.512. The downwash for M = 0 in figure 11(a) is good although improvement for 
y = 0.8 could be desired. In figure 11(b), M is increased to 0.8 and the distribution 
at y = 0.8 is noticeably worsened. Investigation of a range of E s might reveal a bet- 
ter overall choice. 

The discussion now turns from consideration of fundamental or first-order pressure 
and associated downwash modes to consideration of higher-order pressure and downwash 
modes. A commonly used method for developing higher order pressure modes has been 
to multiply the fundamental mode by powers of (cos 0 - cos 0 C ), as in equation (7). It 
was desired to determine what would be the effect of such multiplication. The choice 
here is to use as the reference not the hinge line 0=0 c, but rather 6 = 6 C q for the line 
on the flap midway between its hinge and trailing edge. The very last terms that contain 
C 00 (M) equations (32a) and (33) are of the type described. With the definitions 




bos 6 - cos 6 




cO 1 


M 


Coo(M) 

2 


B 2 0 (M) = 


/b\2 / cos 6 _ ~ cos g c o\ _1_ 

W \ M / “ 12 


|Coo(M) 
J 


(36) 


the downwashes WBjq an( l Wj^q are originated and for the rectangular wing with an 

aspect ratio of 4 are shown in figures 12(a) and 12(b) for M = 0 and Es = 1/2. Inclu- 
sion of the quantities p and 1/12 has only quantitative, not qualitative, effects. The 
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downwashes Wggg of figure 4(c), Wg^g, and Wg 2 g are for the same configuration 
and represent the first three terms of a series that could be combined linearly in an 
effort to fit a specified distribution. Hence, they are to be examined for the desired 
orthogonality. The downwash Wg^g of figure 12(a) shows a chordwise slope and passes 
through zero as is appropriate for a chordwise second-order mode. The downwash 
Wg 20 of figure 12(b) shows relatively little of the curvature which would be appropriate 
for a chordwise third-order mode. Furthermore, Wg^g displays a noticeable resem- 
blance to both Wggg and Wg 2 g. Hence the set is not strongly linearly independent 
and is rather far from being orthogonal (with respect to a unit weighting function). The 

f (O) 

resemblance of Wg^g to of figure 6 is even stronger. This resemblance would 

bring a poor conditioning of the downwash matrix jw] of equation (8b). In addition, the 
capability for rapid convergence of the downwash expansion of equation (8a) (with a small 
number of modes) to the desired downwash distribution over the whole wing, might also 
be impaired. Therefore, the higher order pressure terms of the type B'^q( 0,??) and 
B 2 g(0,?]) are judged to be less than adequate. 

In contrast, use of the finally adopted pressure terms B^g(0,??) from equa- 
tion (32a) and B 20 (d>V) from equation (33) produce the downwash distributions Wgjg 
and Wg 20 of figures 13(a) and 13(b), respectively, for the same wing parameters as in 
figure 12. The distributions on the control, featuring the steep slope and passage through 
zero of Wg^g and the strong curvature of Wggg, are judged to represent a marked 
improvement toward achieving rapid convergence of the downwash series in a collocation 
procedure. The set Wggg of figure 4(c) , Wgjg (fig. 13(a)), and Wggg (fig. 13(b)) 
are not perfectly orthogonal over the flap surface but are probably close enough for the 
downwash fitting process. 

Figures 14(a) and 14(b) show the results of using the pressure modes Bg2(0,f?) 
and Bg4(0,?]) to produce the downwashes Wggg and Wgg^, respectively. The latter 
are primarily spanwise varying second- and third-order modes. In figure 14(a), Wggg 
shows a nearly linear spanwise variation on the control and nearly step -function discon- 
tinuities at the control ends. In figure 14(b), Wgg^ shows a parabolic curvature on the 
control and nearly step-function discontinuities at its ends. The variations of Wggg 
and Wgg 4 are markedly less on the wing than on the control. These characteristics 
are in accord with desired goals. 

CONCLUDING REMARKS 

For partial-span control surfaces of the trailing-edge -flap type with leading-edge 
hinge and sealed gaps, a series of lifting pressure distributions have been developed for 
deflected and deformed flaps on thin, swept, tapered wings for use in aerodynamic and 
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aeroelastic analysis. The development is based on the lifting-surface integral equation 
that relates lift and downwash. An important advance is made by hypothesizing an anal- 
ogy with the flow characteristics over a nonlifting wing with symmetric thickness dis- 
tribution and then applying the determined pressure perturbation characteristics to the 
lifting pressure on a wing with control. 

For lifting wings, the fundamental or first-order pressure term gives a downwash 
distribution that is a good approximation to the following: a small uniform downwash on 
the wing and a step-function jump to a large uniform downwash on the control. A chord- 
wise second-order term gives a downwash with strong linear chordwise variation, and a 
chordwise third-order term gives a downwash with strong chordwise curvature on the 
control; both downwashes are well behaved on the wing. Two higher order spanwise 
varying terms are also developed. These feature downwashes which have strong linear 
spanwise variation and strong spanwise curvature, respectively, on the flap and are well 
behaved on the wing. This type of lifting-pressure series is expected to lead to rapid 
convergence of the downwash collocation procedure for deflected and deformed controls 
and to provide an advanced basis for extension to the oscillating case of flutter. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., May 29, 1969. 
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APPENDIX A 


EXPRESSIONS FOR THE POLYNOMIALS Q r _i(y-i?) 
AND THE FUNCTIONS i/ r (f) 

Expressions for the polynomials Q r _i(y-r)) and the functions v T {\ 
tion (17) for r up to 5 are as follows: 

Q 0 (y-® - jjg 

Qi(y-Tj) = -n) * (y c0 - s) 


r* 2 


Q2(y-h) - 




!(y - 77 ) + |[56 c 2 (y c0 - y ) 2 - l] 


m* - D 2 + ^ 

3/3 2 L 

2 /rrs 2 


Qsiy-v) = 


^ -**♦§ U 6 c 2 (y c0 - y ) 2 


2/32 l 16 


^ - © 2 + ^ 2i2 ] - <y - ® + (y c o ■ y){ 75 c 2 (y c0 - y) 5 

[(x - l) 2 + /3 2 z 2 ] - 3^ 


145, 


3/3" 


Q4(y-^) = 


«c/»*c 4 ,- -,4 «5 3 - 3 26 c 

■(y - vr + — -r— (y - 17) + 


8/3 2 \ 5 


S 2 

3 


i 5 [ 9B c 2 (y c0 - y) 2 - i) 


185 


96 


5/3 2 
2 


M(x - 1) 2 + ^ 2 z 2 ]/(y - v ) 2 + i 5 6 c 2 (y c0 - y)| 3 ^c 2 (y c o - ff 


8/3 4 
206 c 2 


•[(x - 1) 2 + /3 2 z 2 ] - 1/ (y - v) + 36 c 2 (y c o ' y) 2 [i 5 6c 2 (y c0 ~ 


3§ c 2 (y c o - y) 2 - ^ (x - i) 2 + /3 2 z 2 ] + 24 ^ c 


5/34 


-[(x - i ) 2 + /3 2 


from equa- 


y) 2 - io 
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v Q (l) = 1 

= s c(y c0 - y) 

,, 2® = |[ 3B c 2 (y c o - y) 2 - *] - - 1) 2 * 0 2£2 ] 

»,© - %%sc 2 (y c „ - y f - 3] - - »* * ^ 2 ] 

-«® - - foBe^o - y) 2 - 1] 
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APPENDIX B 


EXPRESSIONS FOR THE FUNCTIONS C Hm (x,y, 0 ) FOR 
THE CONFIGURATION IN SKETCH 5 

Expressions for the functions CH 0 Q'(x,y, 0 ), CH 10 (x,y, 0 ), Cji 2 o(x,y, 0 ), 
CHoi(x,y, 0 ), and Cj[Q 2 (x,y, 0 ) for the configuration in sketch 5 are as follows: 


c Hnn(*>y>°) = * ln 


CHio^,y,0) -^(y -y cl ) In 


(y ' ^ 2 ) 


+ /3 2 

(y - y c i) 2 - £ 

(y - y c i) 


+ /3 2 

(y - y C 2) 2 - 13 

(y - y C 2) 



x c ) 2 + A 2 (y - 

yci) 2 + (* • 

x c ) 


A 




Xc) 2 + /3 2 (y - 

yc2) 2+ ( J - 

- x c ) 

A J 


+ =(*c 0 ' x)c Hoo ( 5 ,y, 0 ) 


C H 20 (x,y, 0 ) = -^<(y - y cl )^(x - x c ) 2 + £ 2 (y - y cl f - (y - y c2 )p - x c f + p 2 (y - y c2 ) 


z\ 2 « 2 /r r 


+ 2 (*co - *) (y - y c i) ln 


/(x - x c ) 2 + / 3 2 (y- y cl ) 2 + (x - x c ) 
a 


(y - y C 2) 


/(* - x c ) 2 + £ 2 (y - y c2 ) 2 + (* - x c ) 


+ fe(x c0 - 5 ) 2 - flcH 00 (x,y,o) 

M 


CHoi^ 0 ) = "f ^(* - *c) + 0 2 (y - y c2 J - |(x - x c Y + £ 2 (y - y cl )‘ 


+ B c(y c o - y) c H 00 ( 5 5 y> 0 ) 
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c Ho2 (5 >y> 0 ) 
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APPENDIX G 


NUMERICAL INTEGRATION 

The wing modes of pressure A nm (0,r?) can be integrated by some of the already 
developed and well-known procedures. (See refs. 1 to 6.) For the present calculations, 
the method of integration described in reference 1 is used. This method can also be 
applied to integrate the flap modes of pressure B nm (0,t?) if some modifications are 
introduced into the chordwise integration. 

Each mode B nm (0,Tj) is related to its associated flap mode of downwash 
w B nm ( x >y) through the following expression that results from substituting equations (3) 
and (8a) into equation (1): 


W Bnm (x ’ y) = 


--a? d77 r C nr (0,r7)KrM,x-4(0),s(y-rj)]sin 0 d0 (Cl) 

-1 (y - vr Jo u J 


where the relation between p, m, and r is given by equation (31) and that between 
Bnm ar >d C nr by equation (30). 

Special care is needed in the numerical evaluation of the chordwise integral 



C nr (6,rf)K sin 0 d0 


(C2) 


The main difficulty comes from the logarithmic singularity that the function Coc )(8,rj) 
has along the control- surface hinge. This singularity also appears, through the influence 
of the same Cqo(0> v), in all higher order functions C nr (6,rj). 


In cases of a full- span-flap configuration, the elemental logarithmic singularity 
ln^| 0 - 0 e | sin 9q) can be extracted from the functions C nr (8,rj) and the more simple 
integral 

7T 

f ln(|e-8c| sin 0 C ^K sin 0 d0 | 


can be evaluated in three steps: 


or 



(0 X < 0 C ) S (C3) 
(0 X > 0 C ) 
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where 9 X and 6 C represent the chordwise coordinates of the control point and the 
control-surface hinge, respectively, at the station 77 . The integrals from 0 to 9 X and 
from 6 X to 7 t can be evaluated by a Legendre-Gaussian quadrature, with unity as the 
weighting function. The rest of these integrals, in which 6 C appears as a limit of inte- 
gration, can be evaluated by the Berthod-Zaborowski quadrature (see table 9 in ref. 27), 
which provides highly accurate numerical results. 

But this procedure of integration cannot be applied to the more general and common 
case in which the control- surface spans are smaller than the wing span. Further, it is 
not feasible to divide the wing span in several regions isolating those where the control 
surfaces are located because, in spite of the fact that the function Cqq(9,ti ) is not sin- 
gular at any point outside a control-surface hinge, it has numerical values that increase 
without limit as the hinge end is approached. This effect introduces an important diffi- 
culty in the numerical integration of the spanwise regions which do not have a control 
surface. 


The adopted solution has been to extract (subtract and add) the function 
In Nq( 9, 6 c ,r]) where 


No {e,0c>v) = 


y(cos 9 - cos 6 C ) 2 + j 6 2 s e 2 l 

(y Q i - M 

) 2 + / 3 s e (y c i - hi) 

^(cos 9 - cos 0 C J 2 + /3 2 s e 2 | 

2 

(y.c2 “hl) + ^ s e(y C 2 " M) 


(C4) 


from the expressions C nr (9,ij) and evaluate its chordwise integration in closed form. 
Since the function In Nq^^cjT) is an elementary representation of the singular part of 
the functions C nr ( 0 ,? 7 ), the expressions which remain as a result of this extraction can 
each be integrated numerically without difficulty when suitably arranged. 

The singular part of each function C nr ( 0 ,? 7 ) is defined as In N nr (9,9 c ,ri) where 
for C 00 (e,T 7 ), 


Noo(Mc>^) s n 0 (?7)n(e,e c ,77) 


(C5) 


for C io(9, 7]), 


for C 01 (e, 77 ), 


'k COS 0-COS 0 C Q 


®lo{ e > e c>v) - [ n o( ? i) n (Mc>7?jp bc 




N 0l{ e > e C’V) ® ^ 0 m(9,9 C)V )] 165c ^ c0 ~^ 


(C 6 ) 


(C7) 


and so on. 
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For general nr, from equation (C2), 

FnrOr,*,y) 3 




+ m c^ 


5jr(77j'x,y) 


(C8) 


Then, with the definitions in equations (C5) to (C7), the problem is reduced to the numeri- 
cal integration of the expression 


Fj r (r?;x,y) = j In NnrfMo^K sin 0 d0 
By introducing the function In this expression becomes 

F* r (rj;x,y) = j* H nr (0,?])d0 + J In N o (0,0 c ,?7)k sin 0 d0 


(C9) 


(CIO) 


where 


H (Ov) - In MM j? . „ 

H n r(M)-ln 1 ^^Ksin0 

and where H nr (0,??) is not singular at 0 = 0c* 

With the definition 

<P(0,V) = |^K|M,x-^(0),s(y-?7)J - K[M,x-| c (0 c ),s(y-7?)]J In No(0,0c,7)sin 9 (Cll) 
which is zero at 0 = 0 C , then 

Fn,r(??;x,y) = ^ H n r(0,V)d9 + <K9,rj)d$ + KjM,x-^ c (0 c ),s(y-7 7 )]x(0 c ,7 ? ) (C12) 


where 


^c>v) = 


in ^(cos 0 - cos 0 C ) 2 + /3 2 Se 2 (y c i - \v\f + gs e (y cl ~ N) Bl 
Q ^(cos 0 - cos 0 C ) 2 + ^ 2 s e 2 (y c2 - | rT|) 2 + /3s e (y c2 - |??|) 


sin 0 d0 


(Equation continued on next page) 
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e c ) In 


(1 + cos e c f + /3 2 s e 2 (y c i - \ v\Y + ffs e (y cl - H) 


'(l + cos 0 C ) 2 + /3 2 s e 2 (y c2 - hi) 2 + 0s e (y c2 “ I ^l) 


'(1 - cos 0 C ) 2 + g 2 s e 2 (y c i - hi) 2 + |3s e (y c i - hi) 


+ (l - COS 0 C ) ] 


'(1 - cos °c) 2 + £ 2 Se 2 (y c2 - j V |) 2 + ^s e (y c2 - hi) 


0s e (y Gl - hi) In 


- ^e(y c 2 - hl) ln 


(l + cos 0 C ) 2 + /3 2 s e 2 (y cl - 
|(i - cos 0 C ) 2 + £ 2 s e 2 (y cl - 

|(l + cos 0 C ) 2 + |3 2 s e 2 (y c2 - 
(1 - cos 0 C ) 2 + |3 2 s e 2 (y c2 - 


h|f + (l + cos 0 C ) 

hi) 2 - (i - cos e c ) 

hi) 2 + (l + COS 0 C ) 
hi) 2 - (i - cos 0 C ) 


p 7T p 7T 

The integrals \ H nr (0,?7)d0 and \ 0(0,7?)d0 are evaluated by Gaussian quad- 

-'O J 0 

ratures, with unity as a weighting function, once they are divided into the three intervals 
defined in equation (C3). 


For the planform considered, 


0 C = 77 - arc cos 1 - 1 — - jj — \) 


(C14) 


d x = arc cos 


p[s tan A - |(1 - X)] 

1 - y(i - >0 


(y > 0) (Cl 5) 


When the function X(0 C ,?]) is evaluated at the points | *»? | = Y c l or hi = ^c2’ or 
at any point very close to y c j or y c2 , a minor numerical difficulty appears. It comes 
from the term 


0 s e(y c i -hl) ln 


'(l + cos 0 C ) 2 + /3 2 s e 2 (y cl - hi) 2 + (1 + cos 0 C ) 


(l - cos 0c) 2 + 0 2s e 2 (y c l ' hi) 2 ~ i 1 ~ cos e c) 


and the equivalent term for y c2 in equation (Cl 3). 
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The logarithmic functions are indeterminate at the points j? 7 j ~ y c j and 1 77J = y^, 
respectively, and nearby can produce an ’’overflow" or an "indefinite" number, depending 
on the computer, that causes computation to stop. The adopted preventive measure is to 
specify two narrow bands: 

0 

and ) (Cl 6) 

(yc2-M)- e ^ 

and define as zero therein the value of the products (y c ^ - j 77 |j and (y c 2 - \rj\j multi- 
plied by their respective logarithmic functions. The value of e used in the calcula- 
tions is 10" 

The same difficulty appears in the nonsingular parts of the functions C n o(6,rj) 

(the parts that do not include Cqq) and the same preventive measure is used. 
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Figure L- Continued. 
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Figure 5.- Three-dimensional haff-span view of the downwash distribution W bqq over same rectangular wing as in figure 4 for E s = 1/2 (as in fig, 4(c)). 







Figure 7.- Concluded. 
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Figure 9 .- Continued. 
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Figure 11.- Concluded. 
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